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We investigate finite temperature corrections to the Landauer formula due to electron-electron 
interaction within the quantum point contact. When the Fermi level is close to the barrier height, 
the interaction is strongly enhanced due to semiclassical slowing of the electrons. To describe 
electron transport we formulate and solve a nonlocal kinetic equation for the density matrix of 
electrons. The correction to the conductance G is negative and strongly enhanced in the region 
0.522- < G < 1.02|-. Our results for conductance agree with the so-called "0.7 structure" observed 
in experiments. 
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Introduction. The conductance of a quantum point 
contact (QPC) - a ID constriction in a 2D electron gas - 
has been known to be quantized in units of Gq = 2e^/h 
since 1988 The observed conductance plateaus can 

be easily understood in the single-electron picture 0, . 
Below we write conductance in units of Gq. 

The "0.7 structure" appears on the lowest conductance 
step as a narrow extra plateau at G « 0.7. The structure 
was first observed by Thomas et al. in 19960 and has 
been the subject of numerous experiments since 0,0 S 
The position of the structure varies from 0.4 to 
0.7 depending on the device, and the structure becomes 
more pronounced as the temperature is increased up to 
~ 2K, where thermal smearing becomes significant. It is 
not clear from experiment whether the structure survives 
at T = 0. In a longitudinal magnetic field, which breaks 
the spin degeneracy, the structure evolves smoothly to 
the 0.5 plateau that is expected in the single-electron 
picture. This demonstrates that the effect is in some 
way related to electron spin. 

Observations of the 0.7 structure have created much 
theoretical interest. There have been suggested expla- 
nations based on spontaneous magnetization within the 
contact 0, 0, Il4 Hi! ll^ . charge density waves 
within the contact |l8l Il9j. the Kondo effect [i^, l2l| . 
and even on electron-phonon scattering in the contact 
[2^ . The problem has also stimulated the development of 
general scattering theory in the presence of leads |2^[3|- 

In the present paper, using perturbation theory, we 
consider the finite temperature correction to the Lan- 
dauer formula for the conductance of a QPC. Our results 
are very similar to experimental data on the 0.7 struc- 
ture. 

Model. Near the center of the QPC the single-particle 
dynamics are described by the parabolic saddle-point po- 
tential U = Uq — ^muj'^x^ + ^mujyy^, see, e.g., Ref. 
Here m is the effective mass, C/q is an energy that can 
be controlled in experiment using the gate voltage, and 
LOy and ljj are parameters with typical values in experi- 



ment of bJy/uJ 



4meV. Throughout the paper 



we set fi — ks — ^. The dynamics in the x-direction are 
described by different channels corresponding to quanti- 
zation in the y-direction. The lowest, n = 0, channel has 
longitudinal potential —^muj^x^, and hence the trans- 
mission probability is [2 
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i, = (l + e"2-)-^ , 

where e = E/lo, E is the energy. At e = 
point for the nth channel is xt = ±^ 
barrier, from X2d ~ 2/ ^Jmw, many transverse channels 
are occupied and hence the electron-electron Coulomb 
interaction is strongly screened. The interaction is un- 
screened only on the top of the barrier where the electron 
density is low. Therefore the effective electron-electron 
interaction in the n = channel can be approximated as 



k\xi — X2\ 



H,„t = W7r^ge<5(6)'^(6) : 



(2) 



where e is the electron charge, k is the dielectric constant, 
and f — \Jmiox is the dimensionless distance. In GaAs, 
K « 13, m « 0.07me, and for a; ~ ImeV the dimension- 
less coupling constant g^. ~ -^^^ raj u: is about unity, 
5e 1- Another issue related to the many-body screen- 
ing is the absence of pronounced structures in higher con- 
ductance steps. For higher steps there are always lower 
channels penetrating the QPC. This leads to high elec- 
tron density in the contact and hence to screening of the 
effective interaction (2). 

To represent single-particle wave functions in the n = 
channel we consider a ID wire of length L with a potential 
barrier in the middle of the wire. The details of the po- 
tential shape are unimportant apart from the parabolic 
top, see, e. g., 0|. We need the "wire" to normahze 
wave functions according to \il)k^dx = 1, since this 
normalization is convenient for discussion of the nonlo- 
cal kinetic equation, see below. The length L, which is 



2 



of the order of the balhstic mean free path, disappears 
from final answers. Away from the barrier the wave func- 
tion "i/ifc is the standard combination of incident, reflected 
and transmitted waves. Near the potential top the wave 
function is of the form, see, e. g., Ref. I25I 
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where vp is the Fermi velocity far from the barrier, D,^ is 
the parabolic cylinder function, and i' — ie — ^. Eq. Q 
corresponds to the wave incident from the left, A: > 0. 
Due to semiclassical slowing, the probability density at 
the top of the potential. 



7rexp(7re/2) 



^/2cosh(7^e) |r(3/4 - ie/2)|^ 
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is peaked at e « 0.2, see Fig. Q] This results in enhance- 
ment of the interaction (01, and in the end this leads to 
all effects considered in this paper. 



(0. We have checked numerically that the same is true 
with up to higher orders in perturbation theory. Thus the 
interaction (O treated perturbatively does not influence 
the Landauer formula and does not materially change the 
shape of the transmission coefficient at zero temperature. 

Nonzero temperature requires a consideration of the 
details of equilibration. Eq. ((SJ implies nonlocal equi- 
libration, as the scattering states with s = -1-1 equili- 
brate due to collisions in the left lead and the states with 
s — ~1 equilibrate due to collisions in the right lead. As 
a consequence, the many-body density matrix is diago- 
nal not in the basis of plane waves or standing waves, 
but in the basis of the scattering states JSJ (see also Ref. 
[2^). Therefore, the occupation numbers n^, the diago- 
nal matrix elements of the density matrix, obey the ki- 
netic equation — — . Here r is the relaxation 
time in the leads. This equation is in principle equally 
valid for finite bias, but in the present paper we consider 
only an infinitesimal bias where nok is given by Eq. (jSJ. 
Taking into account electron-electron interaction in the 
QPC, we obtain a collision term in the kinetic equation 
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FIG. 1: The probability density at ^ — versus the dimen- 
sionless energy t. 



The current in state JSJ is jk = j^^te, and hence the 
total current obeys the Landauer formula 



J = 2e 



Ldk 
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-^Jknok = —t,V 



nok 
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Here n^k is the non-equilibrium occupation number un- 
der a small applied voltage V, nj — nf{e) is the equi- 
librium Fermi-Dirac distribution, s — k/\k\ shows the 
direction of current flow in a particular state, and n'j — 

— « (5(e — /i); /i is the chemical potential in units of 

UJ. 

The nonlocal kinetic equation. At zero temperature 
the interaction ||2Jl in the leading order is described by 
the usual direct and exchange diagrams for the single- 
particle scattering amplitude. The corresponding correc- 
tions renormalize the frequency lo in the single particle 
potential and even make the effective potential slightly 
nonparabolic (flattened). Nevertheless this does not ma- 
terially change the profile of the transmission coefficient 
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Here Mkk^k^k^ is the matrix element of the interaction 
This matrix element corresponds to the real transi- 
tion between quantum states, \k,ki) 1^27^3)- There- 
fore in the matrix clement the initial and final wave func- 
tions are given by l(2Jl. This differs from the rule for the 
scattering amplitude, where initial wave functions are 
given by (|3Jl, while final wave functions ar e •i/ 'l \x) — 

(5) 

i'-ki^) (t^'2 Sommerfeld rule, see, e.g., Ref.|23|). 



Using occupation numbers from (O and expanding the 
collision integral up to the first power in the bias V ^ we 
find the integral at small temperature, T <^ uj, 



T'^eVL^ 
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All legs in the matrix element are taken at the Fermi 
surface, so we need only perform summations over the 
directions si, 52,53. Calculating the matrix elements of 



3 



interaction (01 with the wave functions Q we find 

St(rtfe) = — sS{Ek-ujfi) 

ovt. 
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where is given by |0J. Eq. JSJ) leads to the following 
steady-state solution of the kinetic equation ^ 



eV 

Uk ^ Uf + —sS{Ek - uj^) 
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Substitution of this expression into Eq. (jsj , where riQk 
Uk, gives an altered conductance in units of Go 




FIG. 2: Conductance in units of 2e^ /h versus fi (chemical 
potential in units of ui) for different temperatures: the weak 
coupling limit, <^ 1. The uppermost curve corresponds to 
geT = 0, while the lowest is g^T « 0.3K, assuming ^^^^ ~ 1. 






FIG. 3: The 
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correction to the matrix element. 
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Eq. pU|l is justified only if the temperature-dependent 
correction term is small compared to unity. The term is 
due to the current of correlated electrons. It will be most 
significant within the region 0.5^ < G < 1.0^, since 
above G = 0.5 electrons can flow without tunnelling. The 
length L is of the order of the mean free path in the leads, 
so it is most natural to assume that the factor in l|10() 
is of the order of unity, though we cannot exclude some 
dependence of the factor on temperature. In the latter 
case the dependence of the correlated current will be 
modified. 

Note that the correction we have discussed is due to the 
interaction between electrons with opposite spins. The 
interaction between electrons with parallel spins vanishes 
because the exchange diagram exactly cancels out the 
direct one for the contact Hamiltonian (0). Therefore 
under an applied longitudinal magnetic field B we should 
take p^ifi) p2(p)p^(//')- Here ^' = - {2gsfiB)B/uj, 
where gs is the gyromagnetic ratio and /is is the Bohr 
magneton. Since p(fi) is a peaked function, see Fig. ^ a 
magnetic field B ^ uj/{2gsfJ,B) ~ 5T effectively switches 
off the interaction. This gives an explanation for the 
smooth evolution of the 0.7 structure to the 0.5 plateau 
of non-interacting electrons under a magnetic field. 

A set of plots of for different temperatures is shown 
in Fig. 121 Though the result looks quite sensible it is 
obtained for ge ^ 1. However, according to our estimate 
the constant is not small, g^ ^ I, and hence virtual 
rescattering must be taken into account. 

Renormalization of the coupling constant. The lead- 
ing correction to the Born scattering amplitude is given 
by the diagrams shown in Fig. |31 This correction is 



equivalent to renormalization of the coupling constant, 
9e fje + 5ge, whcrc Sg^ip) = 2glK{p) and 

OO OO 
— OO— OO 



£1 + £2 - 2/i 



ei - £2 



Here 6{y) is the step function. The first term in K{ijl) 
(diagram Fig.l^Ji) is logarithmically divergent at ei, £2 — > 
+00. The divergence is a result of the contact approxima- 
tion |5J|. When the energy is large, £ ^ 1, the wave length 
is smaller than the size of the barrier, A ^ 1/ yjmw, and 
the contact approximation fails. To fix the problem, we 
introduce an ultraviolet cutoff A, £1 + £2 < A. Depen- 
dence on the cutoff is weak and we will present all results 
for A = 2. The integrals in (|11|) cannot be calculated an- 
alytically. However numerical integration is very simple 
and we present a plot of K{ijl) in Fig. 21 

Since the coupling constant ~ 1, the second or- 
der correction alone is not sufficient. However in this 
regime the Brueckner approximation '295 usually works 
well. Since the kernel K(fi) is independent of external 
momenta, the Brueckner approximation is equivalent to 
the summation of a geometrical progression, and hence 
the renormalized coupling constant gn is 
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glR{p) , Rifi) = 
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[1 - 2g,K{p)Y 
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Plots of i?(/i) for ge — 1 and ge — 2 are presented in 
Fig. 01 Fig. shows a set of plots of conductance tfj. for 
different temperatures, using Eq. H10|l with g^ — > gnige) 
for ge = 1. 
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FIG. 4: Dashed line: the function K{fi) for the second order 
correction. Solid line: Brueckner correction factor R{i-i) for 
Qe — 1. Long dashed line: Brueckner correction factor -R(/i) 
for Qe = 2. 
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FIG. 5: Conductance in units of2e^/h versus /i (chemical 
potential in units of lu) for different temperatures: the inter- 
mediate coupling limit, ge = 1. The uppermost curve cor- 
responds to T — 0, while the lowest is T ^ 0.3K, assuming 

L ~ ^■ 

The results presented in Figs. El and El are very similar 
to the experimental data on the 0.7 structure. Accord- 
ing to our calculation the exact position of the structure 
depends on the coupling constant g,,'. for small Qe it is 
more like a "0.5 structure" and for ge ~ 1 it is a "0.6-0.7 
structure" . 

In conclusion. Within perturbation theory, we have 
considered transport of correlated electrons through a 
quantum point contact. At zero temperature, the ap- 
proach results in the usual Landauer formula and the 
conductance does not show any structures. To describe 
the current at nonzero temperature we have formulated 
a nonlocal kinetic equation for the occupation numbers. 
A current of correlated electrons scales as temperature 
squared at very low temperatures. The corresponding 
correction to conductance is negative and strongly en- 
hanced in the region 0.5^ < G < 1-0^. We beheve 
that these results are directly relevant to the 0.7 conduc- 
tance structure. For the case of weak coupling, the set 
of plots of the conductance for different temperatures is 
shown in Fig. (21 The corresponding plots for intermedi- 
ate coupling, which is probably the regime most relevant 
to experiment, is shown in Fig. [3 Our model is consistent 
with the experimental behavior of the 0.7 structure under 
a magnetic field: a field smoothly "switches off" the ef- 
fective interaction between electrons. Effects considered 



in the present paper have a very simple physical origin: 
the electron wave function at the barrier and hence the 
electron-electron interaction is strongly peaked when the 
transmission coefficient is slightly higher than 0.5. 
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